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List  2f  Symbols 

length  of  first  layer  of  bar 

length  of  second  layer  of  bar 

acoustic  velocity 

Dirac  delta  function  (0  if  t  5^  t) 

normal  strain 

travel  time 

number  of  terms  in  an  infinite  series 
pressure  (compressive  normal  stress) 
peak  pressure 

mass  density  of  bar  material 

time  parameter  used  with  Laplace  transform 

normal  stress 

time  coordinate 

displacement 

area  of  cross  section 

transmission  coefficient 

product  of  transmission  coefficients 
rodulus  of  elasticity 

Heaviside  function  (0  if  t  <  t  and  1  if  t  >  t) 
impulse  (dimensions  of  pressure  times  time) 

Laplace  transform  of  p 
Laplace,  transform  of  u 

space  coordinates  along  length  of  bar  measured  from  left 
end  of  first,  second,  and  third  layers,  respectively. 
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Abstract 


The  multi-layered  cylinder  made  of  elastic 
materials  and  subjected  to  an  arbitrary  end  excitation 
is  solved  using  the  classical  wave  method,  the  Laplace 
transform  method,  and  the  transmissibility  method. 
Complete  numerical  solutions  are  presented  for  a  pressure 
pulse  having  a  rapid  rise  and  an  exponential  decay. 

Curves  for  displacements,  stresses,  and  accelerations 
at  various  sections  of  bars  of  varying  geometries  are 
presented. 
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STRESS  WAVE  PROPAGATION 
IN  LAYERED  ELASTIC  CYLINDERS 

1.  Introduction.  The  problem  of  a  semi-infinite  solid  under  the 
action  of  a  pressure  pulse  that  is  unifonnly  distributed  over  the  surface 
at  any  instant  of  time  is  formally  equivalent  to  the  problem  of  a  rod 
which  suffers  no  lateral  contraction  and  is  acted  upon  by  a  time-varying 
pressure  at  its  end.  In  this  report  three  methods  of  solution  for  this 
problem  are  discussed;  namely,  the  classical  wave  method,  the  Laplace 
transform  method,  and  the  transmissibility  method. 

Bars  with  one  or  more  elastic  layers  are  considered  with  the 
elastic  properties  varying  from  layer  to  layer.  The  geometries  considered 
are  shown  in  Figure  1. 


p(t) 


Figure  1 
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The  solutions  for  the  bars  of  Figure  la  and  lb  are  given  first  using 
the  classical  wave  method.  The  Laplace  transform  method  is  used  to 
find  the  solution  for  the  bar  of  Figure  Ic.  This  solution  is  shown  to 
reduce  to  the  solution  found  using  the  classical  approach  when  a  ■  oo 
(Figure  la)  and  when  the  second  layer  is  rigid  (Figure  lb) . 

The  transmlsslbillty  method  is  explained  in  conjunction  with  the 
rods  of  Figures  la  and  lb  and  is  then  used  to  solve  the  problem  of 
Figure  Id.  That  solution  is  then  compared  with  the  Laplace  transform 
solution  and  is  also  specialized  to  give  the  solution  for  Figure  le. 

The  pressure  pulses  considered  are  the  impulse  of  Figure  2a,  the 
step  function  (or  Heaviside  function)  of  Figure  2b,  and  the  exponentially 
decaying  function  of  Figure  2c  and  2d. 


Figure  2 
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The  first  two  are  used  to  demonstrate  the  classical  wave  approach  and 
to  explain  the  transmiselbility  method.  The  pulse  of  Figure  2c  is 
used  in  demonstrating  the  Laplace  transform  method.  The  pulse  of 
Figuz*e  2d,  representing  a  practical  sh^e  which  can  be  produced  in  the 
laboratory,  is  used  to  demonstrate  the  power  of  the  transmisslbillty 
method.  Thin  report  concludes  with  graphs  representing  displacement, 
acceler-’tlon  anddiress  at  vairious  locations  along  the  rod  for  the  pulse 
of  Figure  2d. 

2.  Basic  Theory.  The  basic  equation  governing  the  response  of  a 
cylinder  to  a  given  end  disturbance  such  as  an  applied  pressure  is  based 
on  the  equation  of  motion  which  follows  from  Newton's  second  law  and  the 
stress-strain  relation  governing  the  behavior  of  the  assumed  material. 
For  an  elastic  material,  the  normal  stress  0  is  related  to  the  normal 
strain  e  throu^  Hooke's  law 

o  -  Ee  (1) 

where  E  is  Young's  modules  of  elasticity.  The  normal  strain  is 
related  to  the  axial  displacement  u  through  the  illation 


e 


3u 

■SSL 


(2) 


where  x  is  measured  along  the  length  of  the  bar.  The  equation  of 
motion  is  easily  derived  if  we  consider  a  free  body  diagram  of  dx 
length  of  the  bar  as  shown  in  Figure  3.  The  bar  has  the  cross  section 
area  A  and  the  mass  density  p  .  The  normal  stress  assumed  in  the 
conventional  manner  to  Increase  in  the  direction  of  increasing  x. 

The  normal  force  is  given  by  the  product  of  the  uniformly  distributed 


normal  stress  and  the  area 


Figure  3 


or 


Applying  Newton's  second  law  of  motion  we  find 


dc 

dx 


dx  A 


da 


P  — 5  “pa 

at 


(3) 

(U) 


where  a  is  the  acceleration.  This  useful  relation  is  valid  for  any 
continuous  elastic  body  or  to  the  material  within  a  given  layer  of  a 
layered  elastic  body. 

Combining  (l),  (2),  and  (3)  we  obtain 


where  the  subscripts  represent  differentiation  with  respect  to  x  or 
t.  The  constant  c  is  the  acoustic  velocity  for  the  material  and  is 


given  by 


(7) 


The  partial  differential  ^nation  (6)  has  the  general  solution 

u  -  f  (t-  2)  +  g  (t  +  |)  (8) 

where  f  and  g  are  arbitrary  functions  of  (t  -  -)  and  (t  +  g), 
respectively.  The  first  of  these  represents  a  progressive  wave 
moving  in  the  positive  x  direction  while  the  second  represents  a 
regressive  wave  moving  in  the  negative  x  direction.  Associated  with 
the  displacemait  in  (8)  are  the  normal  strain 

e  *  -If*  (t  -  |)  +  ig*  (t  +  f  )  (9) 

the  normal  stress 

a  =  E  u^  3  -  i  f'  (t  -  ^  +  S  g'  (t  +  I  )  (10) 

the  velocity 

V  -  u^  -  f’  (t  -  j  )  +  g’  (t  +  ^  )  (11) 

and  the  acceleration 

a  -  u^^  =  f"  (t  -  g  )  +  g"  (t  )  (12) 

where  primes  represent  differentiation  with  respect  to  the  argument 

(t  -  ~  )  or(t+-  ),  as  the  case  may  be.  The  waves  move  with  the  speed 
c  c 

c. 

3 .  Classical  Wave  Method.  The  classical  approach  to  the  problem  is 
centered  around  a  straight  forward  application  of  prescribed  end  con¬ 
ditions  to  the  governing  equation  (8).  Two  elementary  problems  will 
demonstrate  the  method  and  will  be  useful  later  on  as  comparisons  are 
made  between  the  classical  method  and  the  other  methods  considered. 
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Two  pressure  pulses  will  be  considered.  The  pulse  of  Figure  2a 
is  applied  at  x  -  o  and  is  described  mathematically  by 

CT(o,t)  =  -  p^  T  6(t)  (13) 

where  6(t)  is  the  Dirac  delta  function  defined  by 

o  t  <  0 

6(t)  -  { 00  t  =  0  (11^) 

^  o  t  >  o 

00 

J  6(t)  dt  =  1  (15) 

-CO 

The  product  p^  x  is  an  imp’.ase  and  may  be  denoted  by  the  symbol  I. 
The  negative  sign  in  (13)  indicates  a  conpressive  stress,  Ms  may 
rewrite  (13)  as 

CT(o,t)  -  -  I  6(t)  (16) 

The  piQse  of  Figure  2b  when  applied  at  x  «  o  is  described 
mathematically  by 

CT(o,t)  =  -  p^  H(t)  (17) 

where  H(t)  is  the  Heaviside  function  and  is  defined  by 

H(t)  »  I  °  ^  °  (18) 

^1  t  >  o 

More  generally,  we  may  introduce  the  delay  time  t^  airi  Ite 

0  o  <  t  <  t 

6(t  -  t^)  s  “  I  oo 

o  t>  t 


(19) 
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and 


0  o  <  t  < 

H(t  -  t  )  -  H  (t)  -  {  (20) 

0  o  1  t  >  t 

•  0 

Two  bar  geometries  will  be  considered  here:  the  semi-infinite 
bar  of  Figure  la  and  the  bar  of  finite  length  and  right  end  fixed 
shown  in  Figijre  lb. 

(a)  Semi-infinite  Bar.  Our  first  observation  for  the  semi-infinite 
bar  is  that  waves  moving  to  the  right  are  consistent  with  a  pressure 
pulse  applied  at  x=o  but  waves  moving  to  the  left  are  not.  We  con¬ 
clude,  on  physical  grounds,  that 

g(t  +  ^  )  =0  (21) 


and  thtis 

u-f(t-|)  (22) 

The  problem  is  to  find  the  function  f(t  -  f  )  which  is  consistent 

c 

with  the  end  condition  associated  with  the  applied  pressure  pulse.* 
For  the  impulse  of  (16)  we  are  led  to  try  the  general  stress  relation 


a(x,t)  »  -  I  6(t  -  |)  (23) 

which,  at  least,  satisfies  (l6).  In  terms  o^ displacement,  we  use 
(lO)  to  find 

6(  t  -  |  )  (2U) 

or 

u  -+^H(t-i)  (25) 

where 

OD 

H(t  -  t  )  -  /  6(t  -  t  )  dt 
O''  o 

-® 

* 

It  is  assumed  that  all  bars  treated  are  initially  at  rest. 


(26) 


Equation  (25)  is  an  acceptable  solution  for  u  since  f(t  -  7  ) 

c 

in  (22)  is  of  the  correct  form, 

f(t  -  I  )  -  §  H(t  -  2)  (27) 

The  in^julse  pressure  thus  gives  rise  to  a  step  function  type  dis¬ 
placement.  See  Figure  l*a. 

For  the  step-function  pressure  we  have  (17).  A  likely  expression 
for  the  stress  at  any  section  is 

CT(x,t)  -  -  p^  H(t  -  ^  )  (28) 

which  leads  to 

u  (x,t)  =  +  (t  -  ^  )  H(t  -  ^  )  (29) 

where 

jH(t  -  t^  )  dt  *=  (t  -  t^  )  H(t  -  t^  )  (30) 

This  is  a  ramp  function  indicating  that  the  step  pressure  has  associated 
with  it  a  displacement  which  increases  linearly  from  0  at  the  wave 
front  to  an  indefinitely  large  value  behind  the  front.  See  Figure  Itb. 
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(b)  Finite  Bar  with  Fixed  End.  For  the  finite  bar  vd.th  one  end 
(  X  *  a)  fixed  we  must  find  a  displacement  of  the  form  (8)  which 
satisfies  the  end  conditions 

ct(o,  t)  »  -  p(t)  (31) 

and 

u(a,  t)  =  o  (32) 

In  this  case  waves  irrill  be  reflected  back  from  the  fixed  end  and  the 

fvinctions  f (t  -  —  )  and  g(t  +  —  )  are  both  required.  We  must  consider 
c  c 

the  solution  appropriate  to  the  time  interval  of  interest.  For  example, 
vdien  t  <  a/c  ,  only  the  incident  wave  is  present  since  reflection 
can  not  take  place  imtil  an  incident  pulse  reaches  the  fixed  end. 

During  this  interval  we  have 

c(x,t)  -  -  I  6(t  -  I  )  (33) 

u(x,t)  =  +  H(t  -  f  )  (3U) 

for  th^impiolse.  When  t  <~  ,  the  end  condition  (32)  ass\unes 
irportance  and  we  find 

u(x,t)  -  ^  [  H(t  -  2  )  -  H(t  -  ]  (35) 

33. 

When  t  «  —  ,  there  will  be  reflection  from  the  left  end  and  in 

c  ’ 

order  to  satisfy  both  (3l)  and(32)  ,  we  have 
u(x,t)  -  [  H(t  -|  )  -  H(t  _  H(t  -  .^2^)  ] 


(36) 
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lia 

For  t  <  ~  ,  we  have 

-  |2  [H(t  -  I  )  -  H(t  -  )  -  K(t  -  H(t  -  ) 

(37) 

Other  terms  may  be  added  in  order  to  extend  the  interval  beyond 
t  =  iia/c.  For  the  step-function  pressure  pulse,  (37)  become 

u(x.t)  .  ^  [(t  -|)  *(t  )  .  (t  -i5-^  )  H(t  -is-^  ) 

.  (t  )  H(t  -.21^  ) 

*  (t  -  )  H(t  -  )  ] 

(38) 

Three  important  deductions  should  be  made  from  the  above  develop¬ 
ment.  (1)  Reflected  stresses  are  of  the  same  sign  as  Incident  stresses 
when  a  fixed  end  is  reached  but  reflected  stresses  are  opposite  in  sign 
from  incident  stresses  when  a  free  end  is  reached;  (2)  displacements 
reflected  from  a  fixed  end.  are  opposite  in  sign  to  the  incident  dis¬ 
placements  but  are  of  the  same  sign  as  the  incident  displacements  when 
reflected  from  a  free  end;  and  (3)  the  solution  for  multiple  reflections 
is  tedious  using  the  classical  solution,  requiring  a  skillful  adjust¬ 
ment  of  wave  functions  in  order  to  satisfy  the  given  end  ctxiditions. 

1;.  Laplace  Transform  Method.  The  basic  element  in  the  Laplace 
transform  method  is  the  transform  operator  itself.  The  function  f(t), 
for  exaraple^when  multiplied  by  e®*  dt  and  integrated  from  o  to  oo  is 
defined  as  F(s),  the  Laplace  transform  of  f(t).  Thus, 
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t  [  f(t)]  -  F(s)  -  J  f(t)  e  dt  (39) 

o 

If  the  function  to  be  transformed  is  dependent  on  two  variables,  say 
X  and  t,  then 

00 

f  [^f(x,t)J  =  F(x,s)  «  I  f(x,t)  i  dt  (liO) 

o 

The  transform  of  the  first  and  second  time  derivatives  of  a 
function  f(t)  are  found  by  performing  an  integration  by  parts  as 

OO 

£  [f^  (x,t)J  -  J  f ^  (x,t)  e  dt* 

o 

«■  s  F(x,s)  -  f(x,o)  (Ul) 

and 

^tt  F(x,s)  -  s  f(x,o)  -  f^(x,o) 

m 

Having  established  these  relations  we  now  multiply  each  term  (6)  through 
*  s  t 

by  e  dt  and  integrate  from  o  to  go  to  obtain  the  transformed 
equation 


U^(x,s)  =  l^s^  U(x,s)  -  s  u{x,o) 

c 

-  u^(x,o)  ] 

(ii3) 

where 

5  [u(x,t) 

J  •=  U(x,s) 

(Ut) 

If  at  t»o  the  bar  is  at  rest,  we  have 

u(x,o)  ■  u^(x,o)  ■  0 

and 

(il3)  reduces  to 

U^(x,s)  - 

•2^  U(x,s)  -  0 

(U5) 
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The  partial  differential  equation  in  (x,t)  space  is  thus  reduced  to  an 
ordinary  differential  equation  in  (x,s)  space. 

At  a  fixed  end  (x  »  a),  we  have  in  (x,t)  space 

u(a,t)  “  o  (1;6) 

In  (x,s)  space  this  becomes 

U(a,s)  =  0  (U?) 

At  a  free  end  (x  =  o),  we  have  in  (x,t)  space 

a(o,t)»Eu^  (o,t)  =  o  (li8) 

which  in  (x,s)  space  becomes  • 

U^(o,s)  =  o  (1^9) 

At  an  end  (x  «=  o)  with  a  prescribed  pressure  p(t),  we  have  in  {x,t)  space 

a(o,t)  =  -  p(t)  =  Ei^(o,t)  (50) 


and  in  (x,s)  space 


EU. 


,(o,s)  »  -j*  p(t)  e  dt  =  -P(s)  (5l) 


where  P(s)  is  the  transform  of  p(t).  For  the  impulse  pressure,  step 
pressure,  and  the  exponential  pressure  (Figure  2c)  we  have,  respectively. 


P(3)  >=  I  [  6(t)  e  dt  =  I 
o 

CO 

P(s)  =  Pg  j  H{t)  e  dt  -  p^/s 


(52) 

(53) 


P(s)  -  ^  I  (de  -  be  '’^  )  e  dt 


Cd^X 


d  b  --N 

d-*8  ”  b+s  J 


m 


Now  let  us  apply  the  Laplace  transform  method  to  the  two  layered 
bar  of .  Figure  Ic  as  redrawn  in  Figure  5. 


For  this  problem  we  have  the  acoustic  velocities  and  c^  for 

the  two  layers,  the  moduli  of  elasticity  E~  and  ,  the  mass 

densities  and  with  a  length  a  for  the  finite  portion.  The 
prescribed  end  conditions  are 

a(o,t)  ■  -  p(t)  =  Eu^(o,t)  (55) 

u^(a,t)  =  U^Cojt)  (56) 

U2(oo,  t)  =  0  (58) 

Equation  (55)  prescribes  the  end  pressure,  (56)  prescribes  that  the 
displacements  must  be  the  same  at  the  interface,  (57)  that  the  stresses 
must  be  the  same  at  the  interface  and  (58)  that  no  disturbance  be  felt 
at  the  far  end  of  the  bar.  The  coordinates  and  x^  are  measured 
from  the  left  ends  of  the  respective  layers. 

The  governing  equations  in  (x,s)  space  are 

2 

^Ixx  ^*1’®^  ■  ^  ^1  “  o  (59) 

''l 

2 

^2xx  "  ■S  ^2 
°2 


■  o 


(60) 


I 
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Now  let  us  apply  the  Laplace  transform  method  to  the  two  layered 
bar  of  Figure  Ic  as  redrawn  in  Figure  5. 

- a - >1 

p(t)— e>i  t  ^ 

Figure  5 

For  this  problem  we  have  the  acoustic  velocities  c^  and  for 
the  two  layers,  the  moduli  of  elasticity  and  E^  ,  the  mass 
densities  and  with  a  length  a  for  the  finite  portion.  The 

prescribed  end  conditions  are 

a(o,t)  »  -  p(t)  =  Eu^(o,t)  (55) 

u^(a,t)  =  U2{o,t)  (56) 

U2(oo,  t)  =  o  (58) 

Equation  (55)  prescribes  the  end  pressure,  (56)  prescribes  that  the 
displacements  must  be  the  same  at  the  interface,  (57)  that  the  stresses 
must  be  the  same  at  the  interface  and  (58)  that  no  disturbance  be  felt 
at  the  far  end  of  the  bar.  The  coordinates  and  are  measured 

from  the  left  ends  of  the  respective  layers. 

The  governing  equations  in  (x,s)  space  are 

2 

^ixx  ^1  ■  0  (59) 

'^l 

g2 

'^2xx  - 2  ^2  (^2»®^  “  ° 

‘'2 


I 


« 
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and  find 


^2  “ 


P(s)  , 


^1  °2"1 
"-2  ^=1  ! 


2  ‘'1 


ej  P(s)  r  -j 

1  *  A^J  I 


i  -sa/c 

P(s)  !  e  ^  Cl  -  A^2  ) 


1 


Now  substitute  into  (65)  and  (66)  to  obtain 


(Xj^js)  - 


c,  P(s) 


E  s(l  +  A,  „  e 


-  ^12® 


Ug  (X2,s) ,= 


-sa/c-  -sx  /c_ 

P(s)  e  ^  (1  -  A^2)  e  ^  ^ 

-Jsa/c^ 

E^s  (1  +  A^2  ®  ■^  ) 


Using  the  expansion 


-  2^  ( - 1)  » 


we  rewrite  (7U)  and  (75)  in  the  form 


c,  P(s)  .-■.  n  r-  -8(x.^  +  2an)/c, 

1  Q 

-s(2a  +  2an  -  x^)/c3^  -|  (77) 

*12®  J 


-  16  - 


00 


U, (,,,.)  .  Y  *  W=l)/02 

^  0  (78) 


For  the  special  case  of  a  semi-infinite  bar  of  one  layer  ■  pg  =  p, 

-  Eg  -  E,  and  -  Cg  -  c,  and  we  have  A^g  -  o.  Then  (71*)  and  (75) 

reduce  directly  to 


Ui  (x^,s) 


-sx  /c 

c  P(s)  e 


Es 


Ug  (Xg,s)  = 


c  P(s) 


-s(x.+  a)/c 


Es 


(79) 


(80) 


These  must  be  transformed  finally  from  (x,s)  space  back  to  (x,t)  space. 
This  requires  us  to  first  state  the  exact  nature  of  the  pressure  pulse 
whose  Laplace  transform  is  P(s).  If,  for  example,  we  have  the  impulse 
pressure 

then  from  (52) 

and  we  have 


p(t)  -  I  6(t) 

(61) 

F(s)  -  I 

(82) 

_  -sx  /c 

Ui  (x^,s)  -  ^  e 

(83) 

/  X  cl 

^2  ® 

(81*) 

It  is  easily  shown  that  (e"*'®)'/s  is  the  transform  of  H(t  -  k)  and 
Uj^(Xi,s)  and  Ug(x2,s)  are  the  transforms  of  u^(x^,t)  and  Ug(xg,t). 


Thus, 


\  (x^,t)  -  ^  H(t  -  ^) 


(85) 


I 
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..T  X  +  a 

Ug  (X2,t)  -  |i  H(t  -  -  ) 


The  second  of  these  is  quite  clearly  included  in  the  first  and  we  write 


sinply 


u(x,t)  =.  H(t  -  ~  ) 


This  checks  the  result  found  for  the  same  problem  in  (25). 

For  the  special  case  of  the  finite  bar  of  length  a  and  fixed 
right  end  we  have  E2  »  00  and  A^2  “  Equations  (77)  and  (78) 
reduce  in  this  case  to 

-CD. 

c  P(s)  \  r  -s(x  +  2an)/cL  -s(2a  +  2an  -  x  )/c, 

(88) 


U2(x2,s)  -  0 


When  the  pressure  is  a  step  pulse  we  have 


p(t)  ■  Pp  H(t) 


and  from  (53),  P(s)  •=  P^/s.  Thus 


U(x,s) 


CD 

°  Po 


(  -D"  [ 


-s(x  +  2an)/c  -s(2a  +  2an  -  x)/c  J 


where  the  subscript  on  x,E,  and  c  have  now  been  dropped.  The 


transform  of  e  /s^  is  (t  -  k)  H(t  -  k).  Hence,  we  find  the  inverse 
transform  of  (91)  to  be 


00 


-  18  - 


„(x,t) .  142  y  (-1)”  [<‘  -  H 


2an 


)  H(t 


X  +  2an 


) 


.  (t  -  -.  .»)  H(t  -  t  2.a£_i 


-)] 

(92) 


For  the  n  =  o  and  n  »  1  terms  this  expands  to 
P«°  r .  IT  .  Y  .  2a  -  X 


u(x,t)  -  [(t  -  I  )  H(t  -  I  )  -  (t 

X  +2a 


)  H(t  -  ) 


.  (t  -  iL-S  )  H(t  -  ^  (t  -  )  H(t  -  iiSL^  )1 

c  c  c  c  ^ 


(93) 


This  is  the  complete  solution  if  t  <  ^^/c.  It  checks  (38)  which 
was  derived  using  the  classical  wave  method. 

Having  established  that  the  two  layer  solutions  of  (77)  and  (78) 
reduce  to  those  found  earlier  for  the  single  layer,  we  now  convert  the 
displacements  in  (x,s)  space  to  (x,t)  space.  We  consider  particularly 
the  exponential  pressure  pulse  of  Figure  20  having  the  transform  (51*). 

In  (x,s)  space,  we  have 

-b(x  +  2an)/c, 


00 


Z_ 


12 


)"r. 


'*12' 


U2(x2,s) 


Pq  °1  -  *12)  ° 

(d  -  b)  E^s 


-sa/c^ 

"  4 


.s(2a  .  2an  -  x^)/c".  ^ 

•J  td+s 

(91*) 

n  *  2nac2/ci)/c2 


(  -  A,2)“  e 


X  P-  -  -k.] 
Ld+s  b<-8j 


(95) 


» 


-  19 


[ 


I 


At  this  point  it  is  convenient  to  note  a  great  simplification  in 
finding  the  inverse  transforms  of  expressions  such  as  (9li)  and  (95). 
In  general,  the  inverse  transform  of  P(s)/s  is  given  by 

E  [  P(s)/s  ]  '  J  p(t)  dt  (96) 

0 

We  now  write  (91^)  and  (95)  in  (x,t)  space  using 


l'^[  ^  ]  -  J  (de-<“  .  be-  »‘)  dt  .  -e-“  e  e-^* 


(97) 


The  exponential  terms  in  (9U)  and  (95)  give  rise  to  a  translation  of 
the  multiplicative  function.  Thus 


Ui(xi 


p  c  r 


r  .d(t  -  i_L2£2 


) 


2a-f2an-XL 
-d(t  -  ^  •  ) 


-b(t 


''1 

X.  ■¥  2an 


■  ^12  ® 


) 


b(t 


2a+2an-x, 


4e 


-*12  ° 


(96) 


00 


P  Ct 
^0  1 


(-A  )“ 


+e 


Xg  +  2nac2/c^+ac2/c]^ 


x„  +  2naC5/ci^ac2/oi 
^(t  -  ■  )^  (1  -  hi) 

m 


where,  for  simplicity,  the  Heaviside  functions,  similar  to  those  appearing 
in  (92),  have  been  omitted.  Consequently,  any  term  in  (98)  or  (99)  must  be 
taken  as  zero  whenever  the  parenthetical  expression  in  the  exponent  for 
that  term  is  negative. 
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5-  Transmissibilitv  Method.  The  transmissibility  method  discussed 
here  is  based  on  a  direct  physical  interpretation  of  the  results 
found  in  the  previous  section  using  the  purely  mathematical  technique 
of  the  Laplace  transform.  The  object  is  to  be  able  to  come  to  the  same 
results  in  less  time  and  to  have  a  better  appreciation  of  the  physical 
behavior  of  the  system.  We  begin  by  making  several  general  observations. 

First,  in  equations  (77)  and  (78)  we  note  that  every  term  of  the 
series  on  the  right  hand  side  of  the  equation  is  of  the  form 

U(x,s)=^  e"®^B  (100) 

where  the  constant  B  includes  constants  such  as  While 

this  is  specifically  shown  for  the  bar  of  two  layers,  it  is  generally 
true  for  bars  of  any  number  of  layers.  An  increasing  number  of  layers 
only  introduces  more  constants  such  as  quantity  k 

in  the  exponent  of  (100)  will  be  seen  to  depend  on  c,  u,  x,  and  the 
lengths  of  the  layers. 

Second,  U(x,s)  is  the  transform  of  u(x,t).  Its  general  form, 

shown  in  (100),  is  especially  convenient  since  we  can  easily  show  that 

c  B  t 

u(x,t)  =  ^  delayed  by  k  (101) 


For  example,  if 


p(t)  =  I  6(t) 


(102) 


as  in  (81) 
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' 

Fourth,  since  the  displacement  u(x,t)  is  a  function  of  (t  ” 

'  in  the  general  case,  it  follows  that  normal  stress  a(x,t) ,  being 

proportional  to  the  x  derivative  of  u(x,t),  will  be  opposite  in 
sign  to  the  displacement.  Similarly,  velocity  and  acceleration, 
being  related  to  the  time  derivatives  of  u(x,t)  will  have  the  same 
signs  as  the  displacement.  A  positive  change  in  displacement  is  thus 
associated  with  a  negative  change  in  normal  stress  and  positive 
changes  in  velocity  and  acceleration.  The  term  in  (74),  for  example, 

represents  a  negative  change  in  U^(Xj^,t)  but  will  represent  a  positive 
change  in  a^(x^,t). 

With  these  observations  in  mind,  we  now  study  the  results  for  the 
semi-infinite  bar  of  one  layer  and  the  finite  bar  of  one  layer  with 
right  end  fixed  with  the  goal  of  learning  how  to  predict  the  terms  in 
the  series  associated  with  more  complicated  geometries. 

For  the  semi-infinite  bar  of  one  layer  we  have  the  rejsults  (82) 
and  (84)  which  are  the  same  and  may  be  written  in  the  general  case  as 

-S2C 

U(x,s)  =  °  e  ^  (106) 

This  gives  for  the  impulsive  pressure  p(t)  =  I  6(t), 


u(x,t)  =  Y 

(10?) 

a(x,t)  =  —I  6(t  ~  ^) 

(108) 

% 
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c  P 

u(x,t)  =  “^  (t  -  f)  H(t  -  (109) 


a(x,t)  =  “  Pq  H(t  - 


(110) 


In  this  case,  we  have  only  an  incident  wave  advancing  in  the  x  direction. 

For  the  single  layer,  fixed  end,  bar  we  have  in  the  general  case 
from  (77)  with 


-s(x^-2an)  -s(  2a+2an-x)  I 

U(x,s)  =  2^  c  _  ^  c  j 


Now  expand  this  series  to  obtain 


U(x,s)  =  k  -  e  ° 


-sx  -s(  2a-x)  -s(x^•2a)  -s(4a-x) 


-  e  °  +  e 


-s(x+4a)  -s(6a-x)  -s(x+6a)  -s(8a-x) 


+  e 


=  +  e  ° 


+  •  .  .1 


(112) 


The  first  term  represents  the  incident  wave  [compare  with  (106)].  The 
second  represents  the  wave  which  has  traveled  the  length  of  the  bar 
once  and  in  addition  has  traveled  the  distance  a  -  x.  This  is  clearly 
the  wave  which  has  been  reflected  off  the  rigid  and  and  is  moving  to  the 
to  the  left.  The  third  term  represents  the  wave  which  has  been  reflected 
off  the  fixed  end  once  and  off  the  free  end  once  and  is  moving  toward  the 
right . 


Note  that  a  sign  change  occurs  in  the  displacement  upon  being 
bounced  off  the  fixed  end  but  no  sign  change  occurs  when  bounced  off 
the  free  end.  A  wave  which  has  been  bounced  off  the  fixed  end  twice 
and  the  free  end  twice  (x  +  4a)  will  have  suffered  two  sign  changes 
and  if  positive  in  the  incident  phase  will  still  be  positive.  Normal 
stresses  I  on  the  other  hand,  suffers  sign  changes  when  bounced  off  a 
free  end  but  not  when  bounced  off  a  fixed  end.  Table  I  shows  the  values 


of  B  and  k  in  (100)  associated  with  each  phase  as  related  to  a 
sketch  indicating  the  distance  traveled  by  the  pulse. 


Table  I 

We  now  adapt  this  approach  to  the  general  problem  of  the  two 
layered  bar  which  is  semi-infinite  in  length  and  governed  by  equations 
(77)  and  (78).  The  first  term  in  the  bracket  of  (77)  represents 
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We  see  again  ty  comparing  term  3  with  2,  5  with  4,  and  7  with  6,  that 
bouncing  off  the  free  end  (the  left  end)  multiplies  the  signal  by  a 
factor  of  +1.  Bouncing  off  the  interface,  however,  multiplies  the 
signal  by  a  factor  of  (which  is  equal  to  -1  for  a  rigid 

interface) .  A  signal  which  has  bounced  off  the  interface  three  times 
and  off  the  free  end  twice  is  thus  multiplied  by  a  factor  of  . 

Table  II  shows  the  signal  paths,  delay  times  and  coefficients  for  the 
first  layer  of  the  bar  under  discussion. 


the  right.  The  factor  [l  -  represents  the  intensity  factor  for 

a  signal  vriiich  is  carried  across  the  interface  from  layer  1  to  layer  2. 

VJhen  multiplied  by  the  intensity  of  the  signal  in  layer  1  which  reaches 
the  interface,  we  get  the  actual  intensity  of  the  signal  transmitted 
into  layer  2.  The  first  term  in  (114)  represents  the  signal  which  has 
come  directly  into  layer  2  from  the  incident  wave  in  layer  1.  The 
second  term  represents  the  signal  which  has  been  bounced  (within  layer  l) 
off  the  interface  once  and  off  the  free  end  once  and  is  now  transmitted 
through  the  interface  into  layer  2.  Just  as  the  factor  [l  - 
necessary  to  indicate  the  diminution  of  intensity  of  the  signal  in 
passing  through  the  interface,  so  also  is  a  multiplicative  factor  of 
-\2  necessary  for  each  bounce  off  the  interface  before  entering  layer  2 
and  a  factor  of  (+1)  for  each  bounce  off  the  free  end.  Thus,  a 
signal  which  bounced  off  the  interface  three  times,  off  the  free  end 
twice,  and  then  entered  layer  two  would  have  a  total  factor  of  -[l  - 
as  indicated  by  the  fourth  term  of  (114).  Table  III  indicates  the  relevant 
information  for  the  second  layer. 


^ - 

-  a  - >j 

1 

1 

.  -x. 

k  =  a/c^  + 

k  =  %/c^  + 

k  =  5a/c^  +  Xg/c^ 

- - — 

_ ) 

'■ - > 

< - ! - >. 

B  =  (1-Ai,) 

B  =  -  (l-A^,) 

Table  III 
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6.  Bar  with  Three  Layers .  We  now  use  the  physical  understanding 
gained  in  the  previous  section  as  a  means  of  writing  with  minimum 
formal  mathematics  the  solution  for  a  bar  which  is  composed  of  three 
layers,  the  third  of  which  extends  to  infinity  on  the  right  as  shown 
in  Figure  6. 


•-  b 


p(t) 


E-  .  C- 

E  ,  c 

E  ,  c„ 

r— ( 

1 

2*  2 

3,  3 

^^2 

Figure  6 


Associated  with  the  interfaces  we  will  have  need  for  the  constants 
Ax2  as  defined  by  (70)  and  A23  defined  by 


23 


''3 

'3 

I2 

V 

>1  +  — 

_3 

£.3 


(115) 


We  now  shall  attempt  to  predict  the  possible,  paths  followed  by  a 
signal  before  reaching  a  certain  position  in  the  bar  and  then  state 
the  constants  B  and  k  in  (100)  associated  with  the  terms  of  an 
infinite  series  which  correspond  to  those  paths. 
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Consider  first  layer  1.  A  displacement  signal  entering  this 
layer  may  bounce  off  the  first  interface  or  may  be  transmitted  through 
the  first  interface  and  bounce  off  the  second  interface  back  into  the 
second  layer.  Each  time  it  bounces  off  the  first  interface  it  is  multiplied 
t>y  -A^2  moving  to  the  right,  or  by  +A^2  moving ^o  the  left. 

Each  time  it  passes  through  the  first  interface  moving  to  the  right  it  is 
multiplied  by  time  it  bounces  off  the  second  interface  it 

is  multiplied  by  “A2'^*  Each  time  it  bounces  off  the  free  end  it  is 

multiplied  hy  (+1) .  Each  time  it  passes  through  the  first  interface 

moving  to  the  left  it  is  multiplied  by  (l+A^^)  •  note  further  that 
every  signal  upon  reaching  a  given  point  in  the  first  layer  will  have 
traveled  a  total  distance  of  (2ia+20b+x^)  if  moving  to  the  right  at 
the  time  in  question  where  a  and  3  are  integers  and  a  distance 

I 

(2aa+23b-Xj^)  if  moving  to  the  left  where  again  a  and  3  are  positive 

integers  or  zero.  If,  for  example,  o=2  and  3=1  the  path  may  be  as 
shown  in  Figures  7a  and  7b.  The  path  associated  with  a  given  set  (a, 3) 
may  not  be  unique  as  can  be  seen  in  Figure  7c  and  7d.  In  general,  we 
may  say  that  the  number  of  paths  associated  with  a  given  set  (o,3) 
is-  for  the  +x^  and  -x^  terms, 


^  (a+  3  -  1)1 
3!  (a  -  1)! 


(116) 


For  the  example  just  considered,  we  find  N=2,  corresponding  to  the 
two  cases  shown  in  Figure  ?• 
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l[a)  4a  +  2b  + 

(b)  4a  +  2b  - 

(c)  4a  +  2b  +  X 

1 


k —  a  '  b  '  j 

L.  .  ■  I  _L 


(d)  4a  +  2b  -  Xj^  ^ - '> 

Figure  7 


5 


The  coefficient  to  be  assigned  to  any  term  is  dictated  by  the 
signal  path  and  is  given  by  the  product  of  the  factors  associated  with 
the  penetration  of  each  interface  or  the  bouncing  off  each  interface. 

For  example,  the  path  given  in  Figure  7b  requires  the  coefficient  in  (100) 

B  *  ^23 

The  terms  in  parenthesis  represent  in  order  the  penetration  of  the 
first  interface,  the  bouncing  off  the  second  interface,  the  penetration 
of  the  first  interface,  the  bouncing  off  the  free  end,  and  the  bouncing 


Signal  moving  to  the  right 
— _4<c - b  - J 


a  =0,  3^0,  N=1  |a=0,  3=1,  N=1 


B  =  (1-A.  J 


a  =0,  3=3,  N=1 


1-A,.)  A..  A 


a  =1,  3=0,  N=1 


o  =0,  3=2,  N=1 

ChT- 

B  =  (1-A^3)  aI 

Zj 

4 

B  =  -(1-A^2^  ^ 


o  =1,  3=1,  N=2 


(1-A,.)  A 


a=l,  3=1,  N=2 


(1-A^2)  ^23  ■^12 


a  =1,  3=2,  N=3 


B  ^23 


a  =1,  3=2,  N=3 


o  =1,  3=2,  N=3 


B  = 


a  =1,  3=3,  N=4 


B  =  (i-a_;;2)(i-At,)a 


0=1,  3=3,  N=4 


0  =1,  3=3,  N=4 


B=-(1-A-  )(1-A.JA^3A^-  I  B  =  -(1-A^2*^-V4^2l® 


a  =1,  3=3,  N=4  I  I  0=2,  3=0,  N=1 


Table  V 


I 
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a  =2,  0=1,  N=3 


B=(1-A,,)(1-A,,)a,A2 


a  =2,  0=2,  N=6 


r  4 


B=(1-Ai2)  A23  Ai2 


0=2,  0=2,  N=6 


T,-/-.  .2  v2,  .  ,2 

(1-A^2^  ^23 


a  =2,  P=3#  N=10 


o  =2,  0=3,  N=10 


B=-(I-Aij^-Ai2)A|,A^2 
o  =2,  0=3,  N=10 


B=(l-A^^)(l-A^g)A^^j^^ 


a  =2,  0=3,  N=10 


a  ^2,  0=1,  N=3 


B=(1-A3^2H1-At  ^)A23Ai2 


a  =2,  0=1,  N=3 


ID 


a  =2,  0=2,  N=6 
3 


Z) 


B  -(l-A^^)(l-A^^)Ag^A^^  B=-'(l-A^^)(l-Aj^^)j|^A^^ 


a  =2,  0=2,  N=6 


B=-(1-A^^)(1-Aj^^)A^3A^;^ 


o  =2,  0=3,  N=10 


B-(1-A^2^  ^^”^2^  ^23^12 


0=2,  0=3,  N=10 


B=-(l-A^g)^(l-iL,2^A|^^ 


o  =2,  0=2,  N=6 
.  'D~  '  "D 


a  =2,  0=2,  N=6 


B  -(I-aLX^-S.2^43^2 


a  =2,  0=3,  N=10 
-———3 


B-(  1-A3^2^  ^  ^”^12^  ^23^2 


a  =7,  "0=3,  N-=iU' 


cz:^ 


o  =2,  0=3,  N=10 


B  -(I-A-^2^  ^^“^12^*23^12 


B=(l-A^%)(l-A.)g,)A^.,A3^ 


o  =2,  0=3,  N=10 


D 


<■ 

B=(l-A^^)  (1-Aj^2^^23^2 


Table  V  (Continued) 


k -  a - ik— 

b - ». 

®  =0,  3=1,  N=1 

(— - 1 - i - V 

1 

1 — *-^2 

;■ 

Signals  moving 

to  the  left 

B— (1-A^2)  ^23 

a  =0,  p=:2,  N=1 

o  =0,  3=3,  N=1 

0  =1,  3=1,  N=1 

< — ' 

B=+A^2A23(1"-\2^ 

r'  '  ■■' 

^ 

®“"^2^23^^“'^12^ 

. . 

B=+A^2^23^^"'^12^ 

a  =1,  B=2,  N=2 

0=1,  3=2,  N=2 

- ) 

0=2,  3=1,  N=1 

;=> 

C-_  , 

®""^L2^23^^"^12^ 

0  -1,  &=3,  N=3 

^ _ > 

a  =1,  3=3,  N=:^  ' 

< - ^ 

a  =i,  ti=3,  M=3 
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r 
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0=2,  3=3,  N=^ 
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- ; - ) 

w__,.  ....^ 

B=A^2^|^(1-A3_2)(1-A^2^ 

. 

B=A^gA|3(l-A^2^^^~^2^ 

Table  VI 


Values  of  k  for  each  signal  path  are  found  using  the  expi'dssion 


k  =  UtUli  +  SSk  t  i  (la) 

'l  '2  °2 

Each  signal  path  vrith  its  associated  values  of  B  and  k  may  be 
identified  iri.th  a  term  such  as  (100).  All  such  terms  added  together 
will  yield  the  series  which  represents  the  displacement  at  any  point  in 
the  layer  considered  exsjressed  in  terms  of  (x,s)  space.  Equation  (113), 
it  will  be  recalled,  was  such  an  expression  for  the  first  layer  of  the 
two  layered  bar. 

We  translate  the  series  in  (x,s)  space  into  a  series  for  (x,t) 
space  by  performing  the  inverse  Laplace  transform  as  explained  in  the 
previous  section.  The  general  term  (100)  will  transform,  for  example, 
into  the  term  (101).  For  the  impulse  I6(T)  this  becomes  (104). 

The  normal  stress  is  found  by  multiplying  the  first  x  derivative 
of  u(x,t)  hy  the  value  of  E  in  the  layer.  The  acceleration  is  found 
taking  the  second  t  derivative  of  u(x,t). 

6.  Numerical  Example.  As  a  numerical  example  of  practical 
importance,  we  consider  a  pressure  pulse  as  depicted  in  Figure  8. 

It  was  chosen  to  simulate  a  pressure  -  time  function  actually  reproducible 


p(t) 

Po 


Figure  8 


>t 


The  bars  considered  have  the  geometric  properties  shown  in  Figure  9- 
The  length  b  of  the  second  layer  is  chosen  so  that  the  travel  time  T^ 
of  a  wave  through  that  layer  will  be  exactly  1.5  times  the  travel  time 
of  a  wave  through  the  first  layer.  Thus 

b=l. 530^/0^  (125) 


Displacements  and  stresses  as  functions  of  time  were  determined  at 
locations  Xj^=a/3»  x^=a,  and  x^=h.  Accelerations  were  determined 

at  =  a/3  and  =  b/3-  ^ 

The  physical  properties  of  the  materials  in  the  bar  are  assumed 

such  that  ^^3^°3  ” 

To  fix  ideas,  we  now  choose  a  particular  signal  path  associated 
with  a  particular  point  in  the  bar  and  find  the  related  term  in  the  series 
expansion  for  u(x,t).  We  choose  as  an  example  the  point  x^^  =  a/3  in 
the  bar  of  Figure  9c  and  the  path  associated  with  a  =  2  and  0=1 
and  a  signal  moving  to  the  left.  This  gives  a  value  of  k  in  (100)  of 

k  =  (4a-Xj^)/c^  +  Zhjo^  (126) 


Figure  9 


The  transmission  coefficients  in  this  case  are  found  using  (105)  to  be 

(127) 


Ai2  =  0.2 


A23  -  1.0 

From  Table  VI  we  now  find  B  to  be 

B  =  (1-0.04)(1.0)(0.2)  5  0.192 
The  associated  term  in  (x,s)  space  is 


(128) 


U(x  is)  = 
1 


0.1?2  p(s)  -s(4a/c^+2b/c2-Xj^/cj^) 


"1^ 


(129) 


To  facilitate  matters,  we  now  use  (96)  to  find  the  inverse  transform 
of  P(s)/s  where  p(t)  is  given  by  (122).  This  yields 


.  -1  &  ^  n  Ji.  "^i  ^"'^^1^  r  •  575(  )  /t| 

'  *  0.575  J 


(130) 


'?a*n»o 
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The  inverse  transform  of  (129)  is  then 


u{x  ,t)  =  - 
1 


0.192  jTt-k) 


G^(4T^)  + 


1  Q  < _ 

T,  H(t-4T,-k)  r  -0.575(t-4T^-k)/T:nl 

^  ^  0:575-^'-'  ^ 


where 


■  H(t-4T^-k) 


(132) 


We  now  proceed  to  plot  this  curve.  While  the  expression  (I3I)  may 
appear  to  be  tedious  to  plot  it  will  be  found  with  a  little  practice 
that  it  offers  no  great  difficulty.  In  a  like  manner,  we  calculate  every 
other  term  in  the  expansion  for  u(x,t)  and  plot  it.  The  sum  of  all 
terms  is  easily  found  graphically  by  summing  the  ordinates  of  all  the 
curves. 

The  displacement  at  x^~a.l3  for  the  bar  in  Figure  9c  is  shown  in 
Figure  10.  The  curves  for  several  of  the  typical  terms  in  the  series 
expansion  are  shown  as  is  the  curve  representing  the  sum.  Because  every 
term  was  not  considered  in  arriving  at  the  sum,  the  resultant  curve  is 
not  precise  quantitatively  although  it  clearly  describes  the  correct 
trend.  Figures  11,  12,  13,  and  14  show  the  resultant  displacement 
curves  for  all  the  bars  of  Figure  9  at  x^=a/3,  x^=a,  and 

X2=b. 

The  stress  may  be  found  directly  from  the  displacement  aquation 

(131)  to  be 


-41- 


CT(Xj^.t)  =-  0.192  Gj^(4T^)  + 

'  i  ^ 

-0.575(t-4T^-k)/TJ 
+  H(t-4T^-k)  e  i’ 

J 

Curves  for  this  and  several  other  terms  of  the  series  are  plotted  in 
Figure  I5.  Figures  I6,  17,  18,  and  I9  show  the  resultant  stress  curves 
for  all  bars  of  Figure  9  at  Xj^=a/3,  x^=a,  and  x^^b. 

The  acceleration  is  also  found  from  the  displacement  expression 
(131).  It  is 


0.192  c  p 

_ 1,  ..g 

\ 


,Q 

Ht, 


I  1 


G^(4Ti)  - 


-0.575(t-4TT-k)/T^ 
-0.575  H(t-4T^-k)  e  ^  S 

J 


(134) 


This  equation  as  well  as  similar  equations  for  several  other  terms 
appearing  in  the  series  expansion  for  acceleration  are  plotted  in 
Figure  20.  In  Figures  21  and  22  resultant  acceleration 
curves  are  plotted  for  all  the  bars  of  Figure  9  at  x^=a/3;  and 

Vb/3  •  ■  ; 


Figure  12 
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Figure  14 
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Composite  acceleration  curve  for 
all  bars  at 


Figure  22 


